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Abstract 

Transpiration cooling and liquid metals both provide highly effective heat transfer. Using Darcy’ s law in 
porous media, and the inviscid appro xima tion for liquid metals, the local fluid velocity in these flows 
equals the gradient of a potential. The energy equation and flow region are simplified when transformed 
into potential plane coordinates. In these coordinates the present problems are reduced to heat conduc- 
tion solutions which are mapped into the physical geometry. Results are obtained for a porous region with 
simultaneously prescribed surface temperature and heat flux, heat transfer in a two-dimensional porous bed, 
and heat transfer for two liquid metal slot jets impinging on a heated plate. 
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Characteristic dimension; half width of 
jet or coolant slot, p 
Specific heat of fluid at constant pres- 
sure, kcal/kg K 

Volume flow rate per unit depth, nr/s 
Dimensionless jet widths. Fig. 6(a) . 
Thermal conductivity of porous matrix 
kj,j, or liquid metal k, kcal/m s K 
Outward unit normal vector 
Peclet number, F/a 
Static pressure, N/m 2 

Total amount of heat transferred, kcal/s 
Heat flux, kcal/s m 2 

Surface; coolant exit from porous medium; 

half-spacing between jets, m 
Dimensionless temperature, (t - t to )kgi/ 
q a 2b for porous region with free sur- 
face, (t - t OT )k/<^b for jet, 

(t - 0/(1^ r - t„) for porous bed 
Temperature, a 

x and y velocity components, m/s 
Vector velocity, m/s; U = u/(F/2b) 
Dimensionless coordinates, x/2b, y/2b 
Jet stagnation point 
Rectangular coordinates, m 
Thermal diffusivity, K/pCj, m 2 /s 
Permeability of porous medium, m 2 
Fluid viscosity, N s/m 2 
Fluid density, kg/m 3 

Dimensionless potential, <p/F;.O s = <p s /F 
Velocity potential, for porous medium 

<P = (»o ' pJk/uj <P S = (Po " PsWu, 
m 2 /s 

Stream function, m 2 /s; dimensionless 
T = t/F 

Dimensionless gradient, i(d/dX) + j(3/Sy) 


Subscripts 
r : 

s : 

w : 


Reference value 

On surface where coolant exits from 
porous medium 
Wall 


CO 

0 


1, 2 


At reservoir condition 
On surface where coolant enters porous 
medium 

At wall' 1 or' 2 


INTRODUCTION 


In connection with advanced power producing devices, 
the heat transfer engineer must contend with high 


heat fluxes arising from increasing temperatures of 
working fluids to raise thermal efficiency, or from 
high operating temperatures required in devices 
such as fusion reactors. Two means of effective 
codling are transpiration-cooled porous walls and 
liquid metal coolants. A common feature of these 
flows is that locally the fluid velocity is equal 
to the gradient of a velocity potential. This is 
true for Darcy flow in porous media, and for the 
inviscid flow approximation in liquid-metal heat 
transfer analyses. As a consequence the energy 
equation acquires a simplified form when transform- 
ed into a "natural” system of potential and stream- 
line coordinates. For the present problems the 
solutions in the potential plane are reduced to heat 
conduction solutions which are mapped into the 
physical system. 

The use of the potential plane was developed for 
two- and three-dimensional porous cooling in [1) 
and [2], and for liquid metal heat transfer for 
flow across cylinders in [3). The technique is fur- 
ther developed here to obtain: (l) the shape of a 

two-dimensional porous region to provide a desired 
surface temperature for a specified surface heat 
flux, (2) the heat transfer for a two-dimensional 
porous bed with specified boundary temperatures, 
and (3) the temperature distribution on a uniformly 
heated plate cooled by impinging liquid metal jets. 


ANALYSIS 


Figure 1 shows a two-dimensional porous region with 
coolant supplied through a slot. The upper surface 
receives a uniform heat flux and is to be main- 
tained at a specified temperature. The shape of the 
surface is to be found to meet these restrictions. 
The geometry is in dimensionless form in Fig. 4(a) 
and, as will be discussed, maps into potential 
plane coordinates as shown in Fig. 4(b) . Figure 2 
shows a porous region between solid walls Sp and 
Sg with specified temperature distributions. It 
is desired to obtain the total heat transferred to 
the flow. Figure 5 shows the dimensionless form 
and the mapping into the potential plane. The 
third situation analyzed is the impingement heat 
transfer of the parallel liquid metal, slot jets in 
Fig. 3. The dime nsionless form and mapping into 
the potential plane are in Fig. 6. 

Governing Equations 

The solutions obtained here axe for constant prop- 
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Fig. 1 POROUS LEADING EDGE REGION WITH SPECIFIED 
UNIFORM HEAT FLUX AND TEMPERATURE AT SURFACE 


is governed by the continuity and energy equations, 

V • u = 0 (1) 

KV^t - pCpU • Vt = 0 . (2) 

In the porous medium the velocity is governed by 
Darcy' s equation, 

u = - (<c/|a)Vp (3) 

while in the liquid metal flow (assumed irrotation- 
al and invi'scid) the velocity can be written as the 
gradient of a potential 

u = Vq> (4) 

Define a potential <p = (p. - p)(c/u and Eq. (3) 
becomes the same as Eq. (i). Thus the continuity, 
energy and flow equations have the same form for 
both the porous media and jet flows. 

Boundary Conditions 



The thermal boundary conditions are discussed in 
relation to each problem so only the pressure bound- 
ary conditions are given now. The static pressure 
drop as the fluid accelerates from the reservoir to 
the coolant inlet is usually small ccmpared.with the 
pressure drop in the porous medium. Hence, along 
Sq in Figs. 1 and 2 

P = P 0 = Poo = constant x,y on Sq (5) 

The coolant exit face S is also at a specified 
constant pressure 

p = p g = constant x,y on S (6) 

Similarly the inlet and outlet of the jet are at 
constant pressure. . 

Equations and Boundary Conditions in Dimensionless 
Form 


Fig. 2 POROUS BED WITH NONPOROUS SIDE WALLS MAIN- 
TAINED AT SPECIFIED TEMPERATURES 



Let F be the volume flow rate, per unit depth nor- 
mal to the x,y plane, through the porous medium or 
in one of the jets. For the porous medium F is 
found in the analysis in terms of the pressure dif- 
ference p Q - p„, while for the jet it is specified 
as Sbv^ (rig. 3) . A characteristic velocity is 
F/2b as 2b is used as a reference length in 
Figs, 1 and 2. Using dimensionless variables 


Eqs. (1), (2), and (4) become 

V 1 . U = 0 (7) 

$2t - PeU • VT = 0 (8) 

U = V® (9) 

FTom the definitions of q> and 0, conditions (5) 
and (6) become 

0 = 0 X,Y on Sq (ID) 

® = ® s = constant X,Y on S (11) 


Transformation into Potential and Stream Function 
Coordinates 


Fig. 3 PARALLEL SLOT JETS IMPINGING AGAINST A 
UNIFORMLY HEATED .PLATE • . , . ' 

erties. For heat transfer in a porous medium it is 
assumed that the coolant and porous matrix are in 
sufficiently good thermal contact that locally they 
are at a common temperature. With these limita- 
tions the heat transfer in the porous medium or jet 


Conditions (10) and (11) show that the inlet and 
outlet faces of the porous media are at constant 
potential. From Eqs. (7) and (9) 

V • 15 = V • V® = V 2 ® = 0 (12) 

so that in the flow regions the potential satisfies 
Laplace's equation. Orthogonal to the $ lines 
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(a) DIMENSIONLESS COORDINATES WITH BOUNDARY 
CONDITIONS IN TERMS OF A POTENTIAL 
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(b) REGION IN POTENTIAL 
PLANE. 


ions to map into simple rectangles and strips in 
potential-stream function coordinates in Figs. 4(b), 
5(b), and 6(b). This geometric simplification as 
well as simplifications in the differential equa- 
tion and boundary conditions make it convenient to 
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(a) DIMENSIONLESS PHYSICAL PUNE 


Fig. 4 POROUS LEADING-EDGE REGION WITH COOLANT 
SUPPLIED THROUGH SLOT 



(a) REGION AND POTENTIAL BOUNDARY CONDITIONS IN DI- 
MENSIONLESS PHYSICAL PUNE 
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(b) REGION IN POTENTIAL PUNE 


Fig. 5 POROUS BED WITH SPECIFIED WALL TEMH3RATURE 
VARIATIONS 

are streamlines that satisfy 

V 2 ? = 0 C 13 ) 

The impervious boundaries 23 and 56 in Fig. 4(a), 
and Sn and S 2 in Fig. ,5(a) a$p streamlines. 

For the Jet in Fig. 6(a) 12 and 45 are free stream- 
lines, and the flow must_be along the line of sym- 
metry and the plate so 6IB9 is along a streamline. 
The dividing streamline 38 terminates at the stag- 
nation point. These conditions cause the flow reg- 
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Ib) FLOW REGION IN POTENTIAL PUNE 
Fig. 6 IMPINGING LIQUID METAL SLOT JETS 


obtain solutions in the potential plane and then 
conformally map the results into physical coordi- 
nates. 

The pressure boundary conditions are in terms of 
the potential in Eqs. (10) and (ll) . To express 
the energy equation in terms of potential coordi- 
nates, Eq. (9) is used to eliminate U from Eq. (8). 
The resulting equation has the same form as Eq. (16) 
of (1), and in [1] the details are given of the 
transformation into the potential plane. Using 
Eq. (26) of [1] yields the energy equation, 

£t - S^T _ p, 9T = o (14) 

a} 2 a? 2 3$ 


This is the same equation as for convective heat 
transfer in a parallel plate channel with uniform 
flow in the $ direction and with the channel 
width in the T direction. 

From the definition of the stream function the dif- 
ference between the upper and lower streamlines in 
Figs. 4(b), 5(b), and 6(b) is the volume flow rate 
between them. Since T = V/F, the range of ¥ is 
unity. When the geometries in Figs. 4 and 5 are 
mapped into the potential plane the range in 0 is 
0 to 0 S . Thus the determination of 4 S from the 
mapping yields the volume flow rate. 



= f (P 0 - P.) * 


(15) 
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Leading Edge Region with Unknown Surface Shape 

In Fig. I a porous leading edge region is held by 
two insulated supports. The upper surface is sub- 
jected to a uniform heat flux and it is desired to 
maintain this surface at a uniform temperature t g 
set by design considerations. By not having any 
Portion of the surface below t s , the surface is 
not locally over cooled and coolant is not wasted. 
Thus the boundary conditions along S are both 
uniform heat flux and uniform temperature, and the 
shape of S is to be determined. Along the inlet 
face Sq conservation of energy requires that any 
energy conducted out of the porous medium be ac- 
quired by the incoming fluid. Summarizing these 
conditions: 

(16a). 
(16b) 

lyiQ • Vt = pC p (t Q - tJnQ • u x,y on Sq (17) 
In dimensionless form these became 


kj^ng • Vt = q s = constant"! 

> x,y on S 


S s . VT = 1 \ 

f X, 

(T = T s = constant J 


Y on S 


(18a) 
(18b) 

nQ • OT = Be TnQ • U X,Y on S Q (19) 

Since_ S Q is a line of constant potential n = 
and Eq. (19) becomes 0 


• 7T = Be T |TO| 


X, Y on S„ 


(20) 


This problem can be reduced to a heat conduction 
solution which will yield the shape of S. The 
surface S has potential and temperature both ami-, 
form. This suggests trying a solution where T is 
only a function of t>. Then Eq. (14) becomes 

(d Z T/da 2 ) - Be(dT/da) = 0 (21) 

Integrating gives 

(dT/d®) - Re T = Cq (22) 

To obtain Cq apply the boundary condition (20) at 
Sq where 0 = 0. Using VT = (dT/dO ) reduces 
this to 

( ®/ d ®)|o=o = Vo < 23) 

Applying Eq. (23) to Eq. (22) gives Cq = 0. Then 
by integrating Eq. (22) from 0 to ® s , 

-Be(0 s -0) -[(Po-Ps)/ a HK/u) 

T(0) = T s e = T s e (24) 

There still remains to satisfy the uniform heat 
flux along S, condition (18a). Using VT = 

dT/d® and from Eq. (22) SL = Be T s , yields 
the condition for 4>, I S 

= l/(BeT s ) X,Y on S 


“a 


l S (25) 

The boundary conditions for <t> are summarized in 
Fig. 4(a) and from Eq. (12) ® is a solution to 

Laplace 1 s equation. The shape of surface S has 
to be determined such that ® will be constant and 
have a constant normal derivative along S. As 
shown in [4] and [5] this type of heat conduction 
problem can be solved by conformal mapping by using 
an auxiliary potential derivative plane. The freez- 
ing problem [4] depends on the same boundary value 
problem as the. present case. The conformal trans- 
formations are given in [4] for mapping Fig. 4(b) 
into 4(a); to apply these results here requires 
only a change in notation. Results are given in 
Fig. 7. 


Two-Dimensional Heat Transfer in Porous Bed 

When the porous bed in Fig. 5(a) is mapped into 
Fig. 5(b) the problem becomes convective heat trans- 
fer to uniform flow in a parallel-plate channel of 
unit width. Assu m ing the conformal mapping is known 
between Figs. 5(a) and 5(b), the T(x,y) along Sq 
and Sg are transformed to yield the boundary con- 
ditions in the potential plane 


Tq(®) ®,Y on 12 

Tg(®) 4>,Y on 34 


(26a) 

(26b) 


For convective heat transfer in a channel the axial 
conduction term can be neglected unless Ife is less 
than about 10 so for most situations Eq. (14) can 
be approximated as 

(27) 


(a 2 T/3? 2 ) - Pe(dT/d®) = 0 


Without axial conduction the fluid and solid at the 
entrance of the bed are at the reservoir tempera- 
ture. 


T = 0 O,’? on 14 


(28) 


Equation (27) has the same form as the one dimen- 
sional transient conduction equation. The solution 
in the potential plane is found by using results for 
a slab of constant thickness initially at zero tem- 
perature and then heated in a time dependent manner 
as given by conditions (26). From Eq. (1) on 
page 103 of [6] the temperature distribution is 
given by 


*(*.*) = jr S ne‘ (n2/Pe)li2 ® S in ^ 

n=l 

f e (nZ / Pe),lZ> ‘ (Tq(X) - ( -l) n Tg(A) ]dA (29) 

J 0 

Since there is the same volume flow in each AY 
layer, the average exit temperature from the bed is 


T(®_) 


■/v 


Y)dY 


(30) 


The total heat transferred by the bed to the fluid 
is ®tot = P c p F ( t s " O = P c p F ( t w,r “ ^co) T( ^ s ) * 
Inserting T(® S ,Y) from Eq. (29) and integrating 
to obtain T(® s ) yield m 


Qtot 


pCpF( t w ^ r 


_.-L V' 

t J trRe ^ J 


-(n 2 /Be) 


f 


J 

n=l,3,5... 


(n2fl2/Re)A 

e [Tq(X) - (-l^gWldA (31) 


The F is obtained from Eq. (15) by using the ® s 
found in the mapping. 

Heat Transfer for Two Impinging Jets 

As given in [7] the impinging jet in Fig. 6(a) maps 
into the potential region in Fig. 6(b); this repre- 
sents a uniform flow in the 4> direction. With 
axial conduction neglected the heat transfer can be 
computed as convection to two separate flows in 
channels hq and hg wide. The heat losses along 
the free streamlines are neglected. The normals to 
the free streamlines are along the Y direction so 
that 

(dT/dY) = 0 ®,T on 12 and 45 (32) 
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From symmetry at) out streamline 67 

(St/St) = 0 <S,Y on 67 (33) 

The plate 789 is heated with a -uniform flux q w = 

-k Bt/By | y _o or 

(Bt/BY) = -1 X, Y on 789 (34) 

This must be transformed into the potential plane. 
Using BT/BY = (Bt/BY^By/By) 4- (Bt/B<S)(B$/By) and 
the relations along 789 that Y = 0 and Bj/By = 0 
(since v = 0) yield 

-1 = (Bt/By) | y=0 (By/By) | ¥=0 (35) 

From Eq. (9) and the Cauchy- Riemann equations, U = 
Bn/BX = BY/dY. Condition (35) becomes 

(Bt/BY) = -l/u(4>) ®,Y on 789 (36) 

The boundary conditions a£fi^ summarized in Fig. 6(b). 
The velocity U(o) along 789 needed in Eq. (36) is 
given by the conformal mapping solution for imping- 
ing jets in [7]. 

With these boundary conditions the channels of width 
hi and hg have nonuniform heat addition along 78 
and 89-, and the other boundaries are insulated. 

Since Eq. (27) has the same form as the transient 
conduction equation, the condition (36) can be 
thought of as a heat flux variation with time. In 
[6] the transient temperature solution is given for 
suddenly imposed heating at one surface of a slab 
with the other surface insulated. The temperature 
variation at the heated surface for an imposed flux 
of unit magnitude is G(®). By superposition, for 
a variable heat flux. 


w 


-r 

•'ey 


dG j 
— dx 


U(® - x) dx 


(37) 


Differentiating the G(®) given in [6], page 112, 
the temperature distribution along the heated plate 
becomes, 

i»(.) . 



2bq w 

-(Pem 2 h2)'/ X 


1/2 p 

V rPe ) J x 1 / 2 1 U($ - x) I 

-[Pe(nH-l) 2 h|]/x"| 

+ e dx 


(38) 


where Lu = for 0 < X < X D and lu = h- for 
X D — X < oo, and hi and hg are in [77. The map- 
ping relations in [7] are used to transform T w (®) 
to T W (X). Typical temperature distributions are 
in Fig. 8. 


DISCUSSION 

The shapes of a porous cooled region are shown in 
Fig. 7(a). that simultaneously satisfy the conditions 
of uniform heat flux and temperature along the cool- 
ant exit face. The governing parameter involves the 
over all temperature and pressure differences, and 
the surface heat flux. If the allowable surface 
temperature t s is increased, less coolant flow is 
required. The region thickness is increased there- 
by reducing the flow since the pressure difference 
is fixed. If the surface flux q s is increased, a 
greater flow .is required which corresponds to a 
thinner region. The coolant flow is given in 
Fig. 7(b). 



(a) SHAPES OF REGION FOR VARIOUS VALUES 
OF IMPOSED PHYSICAL CONDITIONS. 



<bt COOLANT VOLUME FLOW AS A FUNC- 
TION OF IMPOSED PHYSICAL CONDI- 
TIONS. 


FIG. 7. - POROUS LEADING EDGE REGION. 

Figure 8 shows the temperature distribution along a 
uniformly heated plate for two impinging slot jets 
(from symmetry only half the geometry is shown). 

The surface temperature goes toward infinity as the 
velocity approaches zero in the backflow region at 
x b 0; this Is a consequence of neglecting heat con- 
duction along the streamlines (Bee Fig. 8 of [3]). 



Fig. 8 DIMENSIONLESS TEMPERATURE ALONG HEATED 
WALL FOR IMPINGING PARALLEL LIQUID METAL SLOT 
JETS, S/b = 3. 
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